Universal glassy dynamics at noise-perturbed onset of chaos. A route to ergodicity 

breakdown. 
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The dynamics of iterates at the transition to chaos in one-dimensional unimodal maps is shown 
to exhibit the characteristic elements of the glass transition, e.g. two-step relaxation and aging. 
The properties of the bifurcation gap induced by external noise, including a relationship between 
relaxation time and entropy, are seen to be comparable to those of a supercooled liquid above a glass 
transition temperature. Universal time evolution obtained from the Feigenbaum RG transformation 
is expressed analytically via g-exponentials, and interpreted in terms of nonextensive statistics. 



I 

o 
o 



> 

00 
(N 
l> 
O 
m 
o 



X3 
O 

o 



X 
J3 



While at present the phenomenology of glass formation 
is in a well-documented advanced stage (J the subject 
remains a prevailing and major theoretical challenge in 
condensed matter physics. In experiments and numerical 
simulations the transition of a liquid into a glass mani- 
fests itself as a dramatic dynamical slowing down where 
the characteristic structural relaxation time changes by 
many orders of magnitude in a relatively small space of 
temperatures. Associated to this process, atypical con- 
nections develop between dynamical and thermodynamic 
properties, such as the so-called Adam-Gibbs relation- 
ship between structural relaxation times and configura- 
tional entropy . These poorly understood connections 
pose very intriguing questions that suggest deep- lying, 
hence generic, physical circumstances which may mani- 
fest themselves in completely different classes of systems 
and therefore are capable of leading to novel but universal 
laws. Here we make a case for this premise by exhibiting 
that glassy behavior is in point of fact present in proto- 
typical nonlinear maps close to the onset of chaos. There 
are clear indications that standard phase-space mixing 
is not entirely fulfilled during glass forming dynamics, 
i.e. upon cooling, caged molecules rearrange so slowly 
that they cannot sample configurations in the available 
time allowed by the process 0| . Naturally, the ques- 
tion arises as to whether under conditions of ergodicity 
malfunction, and, as a final point, downright failure, the 
Boltzmann-Gibbs (BG) statistical mechanics is still ca- 
pable of describing stationary states on the point of glass 
formation or those representing the glass itself. The aim 
of this letter is to make evident that the essential ele- 
ments of glassy behavior are all actually present within 
the neighborhood of the chaos transition in simple non- 
linear dissipative maps. By showing this hitherto uniden- 
tified association we put forward a minimal model for 
glass dynamics that endorses the idea of universality in 
this phenomenon, and at the same time provides the rare 
opportunity for detailed examination of properties, such 
as the connection between dynamics and statics and the 
departure from BG statistics. 

A distinctive feature of supercooled liquids on ap- 
proach to glass formation is the development of a two- 
step process of relaxation, as displayed by the time evo- 



lution of correlations e.g. the intermediate scattering 
function i^fe 0- This consists of time t power-law 
decays towards and away from a plateau, the duration 
tx of which diverges also as a power law of the differ- 
ence T — as the temperature T decreases to a critical 
value y 0. This behavior is displayed by molecu- 
lar dynamics simulations Q and successfully reproduced 
by mode coupling (MC) theory Another important 
feature of the dynamic properties of glasses is the loss of 
time translation invariance. This is referred to as aging 
0, and is due to the fact that properties of glasses de- 
pend on the procedure by which they are obtained. The 
slow time decline of relaxation functions and correlations 
display a scaling dependence on the ratio t/t^ where 
is a waiting time. Interestingly, a recent important de- 
velopment is the finding that both two-step relaxation 
and aging are present in fixed-energy Hamiltonian sys- 
tems of N classical XY spins with homogeneous but suf- 
ficiently long-ranged interactions. In these systems the 
length of the plateau diverges with infinite size N oo 
and aging, similar to that found for short-ranged inter- 
action spin glasses, is observed in the long time behavior 
of the autocorrelation function of system trajectories 
. Remarkably, here we demonstrate that the same two 
features, two-step relaxation and aging, are displayed by 
simpler systems with only a few degrees of freedom, such 
as one-dimensional dissipative maps at the edge of chaos 
in the presence of stochastic noise, the amplitude a of 
which plays a role parallel to T — Tg in the supercooled 
liquid or 1/7V in the spin system. This adds to the idea 
of universality for the phenomenon of glass formation. 

The experimentally observed relaxation behavior of 
supercooled liquids is effectively described, under sim- 
ple heat capacity assumptions, by the above-mentioned 
Adam-Gibbs equation, = Aexp{B /TSc), where the 
relaxation time can be identified with the viscosity, 
and the configurational entropy Sc relates to the number 
of minima of the fluid's potential energy surface (and A 
and B are constants) Although at present a first 
principles derivation of this equation is lacking, it offers 
a sensible picture of progressive reduction in the number 
of configurations that the system is capable of sampling 
as T — T„ ^ as the origin of viscous slow-down in 
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supercooled liquids. As a parallel to the Adam-Gibbs 
formula, we show below that our one-dimensional dissi- 
pative map model for glassy dynamics exhibits a relation- 
ship between the plateau duration , and the entropy Sc 
for the state that comprises the largest number of (iter- 
ate positions) bands allowed by the bifurcation gap - the 
noise- induced cutoff in the period-doubling cascade ■ 
This entropy is obtained from the probability of chaotic 
band occupancy at position x. 

So, our purpose here is to contrast the dynamics of a 
fluid near glass formation with that of iterates in simple 
nonlinear maps near the edge of chaos when subject to 
external noise. We chose to illustrate this by consider- 
ing the known behavior of the logistic map under these 
conditions The idea in mind is that the course of 

action that leads to glass formation is one in which the 
system is driven gradually into a nonergodic state by re- 
ducing its ability to pass through phase-space-filling con- 
figurational regions until it is only possible to go across 
a (multi) fractal subset of phase space. This situation 
is emulated in the logistic map with additive external 
noise, Xt+i = ff,{xt) = 1 - nx'^ + 6^, -1 < a;t < 1, 
< ^ < 2, where is Gaussian-distributed with average 
(66') = ^t.t', and cr measures the noise intensity 10] 11]. 
As is well known ^3 CSj the absence of noise a — 
the Feigenbaum attractor at /i = A'c(O) = 1.40115... is 
the accumulation point of both the period doubling and 
the chaotic band splitting sequences of transitions and it 
marks the threshold between periodic and chaotic orbits. 
The locations of period doublings, at /i = < Aic(O), 
and band splittings, at /i = /!„ > /ic(0), obey, for large 
n, the power laws — /ic(0) ^ 6^" and /j,c(0) — /!„ ^ (5"", 
where S = 0.46692... is one of the two Feigenbaum's 
universal constants. The 2nd, a = 2.50290... measures 
the power-law period-doubling spreading of iterate posi- 
tions. All the trajectories with yLtc(O) and initial condi- 
tion — 1 < Xin < 1 fall, after a (power-law) transient, 
into the attractor set of positions with fractal dimen- 
sion df — 0.5338.... Therefore, these trajectories rep- 
resent nonergodic states, as t — > cx) only a Cantor set 
of positions is accessible within the entire phase space 
— 1 < a; < 1. For a > the noise fluctuations smear the 
sharp features of the periodic attractors as these broaden 
into bands similar to those in the chaotic attractors, but 
there is still a sharp transition to chaos at /ic(o') where 
the Lyapunov exponent changes sign. The period dou- 
bling of bands ends at a finite value 2^^°'^ as the edge 
of chaos transition is approached and then decreases in 
reverse fashion at the other side of the transition. The 
broadening of orbits with periods or bands of number 
smaller than 2^^'^' and the removal of orbits of periods 
or bands of number larger than 2^^°"^ in the infinite cas- 
cades introduces a bifurcation gap with scaling features 
[To| [m that we shall use below. When cr > the tra- 
jectories visit sequentially a set of 2" disjoint bands or 
segments leading to a cycle, but the behavior inside each 
band is completely chaotic. These trajectories represent 
ergodic states as the accessible positions have a fractal 



dimension equal to the dimension of phase space. Thus 
the elimination of fluctuations in the limit a —> leads 
to an ergodic to nonergodic transition in the map and we 
contrast its properties with those known for the molec- 
ular arrest occurring in a liquid as T Tg. Since the 
map clearly differs from a Hamiltonian model for a liquid 
in that it does not take into consideration its molecular 
nature we may gain information regarding universality in 
the processes of glass formation. 

The main points in the following analysis are: 1) 
The quasi-stationary trajectories followed by iterates at 
/ic(cr) are obtained via the fixed-point map solution g{x) 
and the first noise perturbation eigenfunction Ga{x) of 
the RG doubling transformation consisting of functional 
composition and rescaling, R/(a;) = af{f{x/a)). Posi- 
tions for time subsequences within these trajectories can 
be expressed analytically in terms of the g-exponential 
function expq (x) = [l — {q—l)x]^^^~'^. 2) The two-step re- 
laxation occurring when tr — > is determined in terms of 
the bifurcation gap properties, in particular, the plateau 
duration is given by the power law txiu) a'^~^ where 
r ~ 0.6332 or r — 1 ~ —0.3668. 3) The map equivalent 
of the Adam-Gibbs law is obtained as a power-law re- 
lation too ~ S-'^, C = (1 - r)/r ~ 0.5792, between t^(cr) 
and the entropy Sc{cr) associated to the noise broadening 
of chaotic bands. 4) The trajectories at //c(o' 0) are 
shown to obey a scaling property, characteristic of aging 
in glassy dynamics, of the form Xt+t^ — h{tw)h{ t/t^) 
where tu, is a waiting time. 

The dynamics of iterates for the logistic map at the 
onset of chaos /ic(0) has recently been analyzed in de- 
tail It was found that the trajectory with initial 
condition Xin = (see Fig. 1) maps out the Feigen- 
baum attractor in such a way that (the absolute val- 
ues of) succeeding (time-shifted t = t + 1) positions ali- 
form subsequences with a common power-law decay of 
the form t''^/'^'1 with g = l-ln2/lnQ;~ 0.24449. That 
is, the entire attractor can be decomposed into position 
subsequences generated by the time subsequences t — 
{2k + 1)2", each obtained by running over n = 0, 1, 2, ... 
for a fixed value of fc = 0, 1, 2, ... Noticeably, the positions 
in these subsequences can be obtained from those belong- 
ing to the 'super-stable' periodic orbits of lengths 2", i.e. 
the 2"-cycles that contain the point x = at /!„ < /ic(0) 
|lOj . Specifically, the positions for the main subsequence 
k = 0, that constitutes the lower bound of the entire tra- 
jectory (see Fig. 1), were identified to be X2" = dn = 



where d„ = 



A2" 



\0) 



is the '71-th diameter' de- 



fined at the 2"-supercycle [lOj. The main subsequence 
can be expressed as the g-exponential xt = exp2_^(— Agi) 
with Xq ~ In a/ In 2, and interestingly this analytical re- 
sult for Xt can be seen to satisfy the dynamical fixed- 
point relation, h{t) = ah{h{t/a)) with a = 2^/'-^-'}^ ,1^ . 
Further, the sensitivity to initial conditions 6 obeys the 
closely related form 6 = expg(A<jt) These properties 
follow from the use of Xin = in the scaling relation 
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that in turn is obtained from the n — > oo convergence of 
the 2"th map composition to {—a)~"g{a"'x) with a = 
21/(1-9). When Xin — one obtains in general 



(2fe 



(2- 



)(0) 



(2) 



When the noise is turned on (cr always small) the 2"th 
map composition converges instead to {—a)~"'[g{a"'x) + 
^aK^G\ (a^x)], where k a constant whose numerically 
determined |14| . value k ~ 6.619 is well approx- 
imated by = 2V2a{l + I/q!^)"^/^, the ratio of the 
intensity of successive subharmonics in the map power 
spectrum |15| . |lClj| . The connection between k and the 
cr-independent u stems from the necessary coincidence 
of two ratios, that of noise levels causing band-merging 
transitions for successive 2" and 2"+^ periods and that 
of spectral peaks at the corresponding parameter values 
fin and Hn+i [isl, |1Q|- Following the same procedure as 
above we see that the orbits Xt at fJ-dcr) satisfy, in place 
of Eq. C^l, the relation 



-1/1-9 



(3) 

where r = 1 — In 2/ In k ~ 0.6332. So that use of Xin = 
yields Xr = t-^^^-i \ 1 -|- (ar^/^^'' I or 



(4) 



where t = t — I and — Inn/ In 2. 

At each noise level a there is a 'crossover' or 'relax- 
ation' time tx — Tx — 1 when the fluctuations start sup- 
pressing the fine structure imprinted by the attractor on 
the orbits with Xm = 0. This time is given by = <J^~^, 
the time when the fluctuation term in the perturbation 
expression for Xr becomes cr-independent and so unre- 
strained, i.e. Xr^ = Tx^^^~'^ \l + i\- Thus, there are 
two regimes for time evolution at ^c{<j)- When t < Tx 
the fluctuations are smaller than the distances between 
adjacent subsequence positions of the noiseless orbit at 
/Xc(0), and the iterate positions in the presence of noise 
fall within small non overlapping bands each around the 
cr = position for that r. In this regime the dynamics 
follows in effect the same subsequence pattern as in the 
noiseless case. When t ^ Tx the width of the fluctuation- 
generated band visited at time Tx = 2^ matches the 
distance between two consecutive diameters, — rfw+i 
where N ~ — Iuct/Iuk, and this signals a cutoff in the 
advance through the position subsequences. At longer 
times T > Tx the orbits are unable to stick to the fine 
period-doubling structure of the attractor. In this 2nd 
regime the iterate follows an increasingly chaotic trajec- 
tory as bands merge progressively. This is the dynamical 
image - observed along the time evolution for the orbits 
of a single state ^ic{<y) - of the static bifurcation gap first 
described in the map space of position x and control pa- 
rameter lU, [l3|, [III. 

In establishing parallels with glassy dynamics in su- 
percooled liquids, it is helpful to define an 'energy land- 



scape' for the map as being composed by an infinite num- 
ber of 'wells' whose equal- valued minima coincide with 
the points of the attractor on the interval [—1,1]. The 
widths of the wells increase as an 'energy parameter' U 
increases and the wells merge by pairs at values Ujq such 
that within the range Un+i < U < Un the landscape is 
composed of a set of 2^ bands of widths Wm{U), m = 
1,...,2^. This 'picture' of an energy landscape resem- 
bles the chaotic band-merging cascade in the well-known 
(x, /J,) bifurcation diagram [l^. The landscape is sampled 
at noise level a by orbits that visit points within the set 
of 2^ bands of widths Wm{U) ~ cr, and, as we have seen, 
this takes place in time in the same way that period dou- 
bling and band merging proceeds in the presence of a bi- 
furcation gap when the control parameter is run through 
the interval < < 2. That is, the trajectories starting 
at Xin = duplicate the number of visited bands at times 
T — 2", n = 1,...,N, the bifurcation gap is reached at 
Tx — 2^ , after which the orbits fall within bands that 
merge by pairs at times r — 2^+", n = 1,...,N. The 
sensitivity to initial conditions grows as = exp^{Xqt) 
(q = 1 — ln2/lna < 1 as above) for t < tx, but for t > tx 
the fluctuations dominate and grows exponentially as 
the trajectory has become chaotic and so one anticipates 
an exponential ft (or q = 1). We interpret this behav- 
ior to be the dynamical system analog of the so-called 
a relaxation in supercooled fluids. The plateau duration 
tx 00 as a 0. Additionally, trajectories with initial 
conditions Xin not belonging to the attractor exhibit an 
initial relaxation stretch towards the plateau as the orbit 
falls into the attractor. This appears as the analog of the 
so-called /3 relaxation in supercooled liquids. 

Next, we proceed to evaluate the entropy of the orbits 
starting at Xi„ = as they enter the bifurcation gap at 
tx{cr) when the maximum number 2''^ of bands allowed 
by the fluctuations is reached. The entropy Sc{fJ,c,tx) 
associated to the state at /Xc(cr) at iteration time tx{cr) 
has the form SdUcjtx) = 2^ as, since at txia) each of 
the 2^ bands contributes with an entropy crs where s = 
— j\pi£,) lnp(^)d^ and where p(^) is the distribution for 
the noise random variable. In terms of tx (since 2^ = 
I + tx and cr = (1 + tx)^^^^^''~) one has Sc^Hctx)/ s — 
(1 + tx)^"^^^^^ or, conversely. 



tx = (S/5,)(1-'-)/'-. 



(5) 



Since tx ~ a''-^, r-1 ~ -0.3668 and (1 - r)/r ~ 0.5792 
then tx 00 and —> as cr ^ 0, i.e. the relaxation 
time diverges as the 'landscape' entropy vanishes. We 
interpret this relationship between tx and the entropy Sc 
to be the dynamical system analog of the Adam-Gibbs 
formula for a supercooled liquid. 

Finally, we draw attention to the aging scaling prop- 
erty of the trajectories xt at ^c(cr)- The case ct = is 
more readily appraised because this property is, essen- 
tially, built into the same position subsequences Xr = 
|c,M(0)|, T ^ {2k + 1)2", k,n = 0,1,... that we have 
been using all along. These subsequences are relevant for 
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the description of trajectories that are at first held at a 
given attractor position for a waiting period of time 
and then released to the normal iterative procedure. We 
chose the holding positions to be any of those along the 
top band shown in Fig. 1 for a waiting time ty^ — 2k + I, 
k = 0,1,.... Notice that, as shown in Fig. 1, for the 
Xin — orbit these positions are visited at odd iteration 
times. The lower-bound positions for these trajectories 
are given by those of the subsequences at times (2fc-|-l)2" 
(see Fig. 1). Writing r as r = + t we have that 
t/t^ = 2" - 1 and xt+t^ = g(*™)(0)5(*/*™)(0) or 

xt+u. =g^'-H0)expg{-Xgt/t^). (6) 

This property is gradually modified when noise is turned 
on. The presence of a bifurcation gap limits its range of 
validity to total times tw + t < (u) and so progressively 
disappears as a is increased. 

In summary, we have shown that the dynamics of 
noise-perturbed logistic maps at the chaos threshold ex- 
hibit the peculiar features of glassy dynamics in super- 
cooled liquids. These are: a two-step relaxation process, 
association between dynamic (relaxation time) and static 
(configurational entropy) properties, and an aging scal- 
ing relation for two-time functions. Our findings clearly 
have a universal (in the RG sense) validity for the class 
of unimodal maps studied. The occurrence of this novel 



analogy may not be completely fortituous since the limit 
of vanishing noise amplitude a —^ (the counterpart of 
the limit T — Tg ^ in the supercooled liquid) entails 
loss of ergodicity. The incidence of these properties in 
such simple dynamical systems, with only a few degrees 
of freedom and no reference to molecular interactions, 
suggests an all-encompassing mechanism underlying the 
dynamics of glass formation. As established 13], the 
dynamics of deterministic unimodal maps at the edge of 
chaos is a bona fide example of the applicability of nonex- 
tensive statistics. Here we have shown that this noner- 
godic state corresponds to the limiting state, a ~^ 0, 
tx oo, for a family of small a noisy states with glassy 
properties, that are conspicuously described for t < tx 
via the g-exponentials of the nonextensive formalism '13| . 
The fact that these features transform into the usual BG 
exponential behavior for t > provides a significant 
opportunity for investigating the crossover from the or- 
dinary BG to the nonextensive statistics in the physical 
circumstance of loss of mixing and ergodic properties. 

Acknowledgments. I am grateful to Piero Tartaglia 
for guiding me into the concepts of glass formation and 
also for his kind hospitality at Dipartamento di Fisica, 
Universita degli Studi di Roma " La Sapienza" . I thank 
Fulvio Baldovin for contributing the figure. Work par- 
tially supported by CONACyT grant P-40530-F. 



[1] For a recent review see, P.G. De Benedetti and F.H. Still- 

inger. Nature 410, 267 (2001). 
[2] P.G. De Benedetti, Metastable Liquids. Concepts and 

Principles (Princeton Univ. Press, Princeton, 1996). 
[3] W. Kob and H.C. Andersen, Phys. Rev. E 51, 4626 

(1995). 

[4] W. Gotze and L. Sjogren, Rep. Prog. Phys. 55, 241 
(1992). 

[5] See, for example, J. P. Bouchaud, L.F. Cugliandolo, J. 
Kurchan and M. Mezard, in Spin Glasses and Random 
Fields, A. P. Young, editor (World Scientific, Singapore, 
1998). 

[6] V. Latora, A. Rapisarda and C. Tsallis, Phys. Rev. E64, 
056134 (2001). 

[7] M.A. Montemurro, F. Tamarit and C. Anteneodo, Phys. 

Rev. E67, 031106 (2003). 
[8] A. Pluchino, V. Latora and A. Rapisarda, 



|cond-mat/03 03081||cond-mat/0306374| 

[9] G. Adam and J.H. Gibbs, J. Chem. Phys. 43, 139 (1965). 
[10] See, for example, H.G. Schuster, Deterministic Chaos. 

An Introduction, 2nd Revised Edition (VCH Publishers, 

Weinheim, 1988). 
[11] J. P. Crutchfield, J.D. Farmer and B.A. Huberman, Phys. 

Rep. 92, 45 (1982). 
[12] C. Beck and F. Schlogl, Thermodynamics of Chaotic Sys- 
tems (Cambridge University Press, UK, 1993). 
[13] F. Baldovin and A. Robledo, Phys. Rev. E 66, 045104(R) 

(2002); cond- mat/0304410 
[14] J. Crutchfield, M. Nauenberg and J. Rudnick, Phys. Rev. 

Lett. 46, 933 (1981) 
[15] B. Shraiman, C.E. Wayne and P.G. Martin, Phys. Rev. 

Lett. 46, 935 (1981). 



X 

X - 

10-' r 



10 




■ •■•.•■•.■■•■•.•■V.VX1WM! 



80, 
48 • ' 



10 



10 



10 



100 



1000 



FIG. 1: Absolute values of positions in logarithmic scales of 
the first 1000 iterations r for a trajectory of the logistic map 
at the onset of chaos //c(0) with initial condition Xin = 0. 
The numbers correspond to iteration times. The power-law 
decay of the time subsequences described in the text can be 
clearly appreciated. 



